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ABSTRACT 

Recent numerical relativity simulations have shown that the emission of grav- 
itational waves at the merger of two black holes gives a recoil kick to the final 
black hole. We follow the orbits of a recoiling supermassive black hole (SMBH) 
in a fixed background potential of a disk galaxy including the effect of dynamical 
friction. If the recoil velocity of the SMBH is smaller than the escape velocity of 
the galaxy, the SMBH moves around in the potential along a complex trajectory 
before it spirals into the galactic center through dynamical friction. We consider 
the accretion of gas onto the SMBH from the surrounding ISM and estimate the 
X-ray luminosity of the SMBH. We find that it can be larger than 3 x 10^^ erg"'^ 
or the typical X-ray luminosity of ultra-luminous X-ray sources, when the SMBH 
passes the galactic disk. In particular, the luminosity could exceed ~ lO^^ergs"^, 
if the SMBH is ejected into the galactic disk. The average luminosity gradually 
increases as the SMBH spirals into the galactic center. We also estimate the prob- 
ability of finding recoiling SMBHs with X-ray luminosities of > 3 x 10^^ erg~^ in 
a disk galaxy. 

Subject headings: black hole physics — ISM: general — galaxies: nuclei — X- 
rays: general 



1. Introduction 



Thanks to recent breakthroughs in numerical relativity, it has been shown that the loss 
of linear momentum radiated away in the form of gravitational waves induces a large recoil 
velocity of the merged binary black hole. This would have been happened for supermassive 
black holes (SMBHs) at the centers of galaxies, if two SMBHs coalesce after a major galaxy 



merger. Since the max imum velocity would reach ~ 4000 km s ^ (e.g. iGonzalez et al.l 12007 



Campanelli et al.ll2007l ). the SMBH could escape from its host galaxy. However, if the recoil 
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velocity is a little smaller than the escape velocity of the galaxy, the SMBH would orbit in 
the potential well of the galaxy for a long time. The identification of observational signatures 
of such recoiling SMBHs is important for studies about the growth of black holes as well as 
the general relativity. 

In addition to direct detection of gravitational waves, a number of ideas have been pro- 
posed for de tection of observatio nal signatures of recoiling black holes through electromag- 
netic waves. iKapoorl (119761 . 119831 ) indicated that stellar captures can lead to the formation of 



an accretion di sk-star system about the SM BH, and t h at the emission from the SMBH could 
be observable. iMadau &: QuataertI (120041 ) and iLoebl (120071 ) argued that a recoiling SMBH 
would be observed as an off- nuclear quasar unt il the gas carried by the SMBH is depleted, 
although iBonning. Shields. &: Salvianderl (120071 ) found no convincing evidence for recoiling 
SMBHs carrying accretion disks in the Sloan Digital Sky Survey data. iMerritt et al.l (120041 ) 
indicated that the displacement of a recoiling SMBH transf ers energy to th e stars in the 
galactic nucleus and converts a steep density cusp into a core. IVolonteril (120071 ) discussed the 
influence of the merger and ejection of SMBHs from the gala ctic centers on the rela t ion o f 
the black hole mass and the velocity dispersion of the galaxy. iGualandris &: MerrittI (l2008l ) 
indicated that helical radio structures could be observed aro und a recoiling SMBH because 
of the oscillation of the SMBH in the core of the host galaxy. iLippai. Frei. fc HaimanI (120081 ) 
showed that prompt shocks are created in the gas disk around a recoiling SMBH and that 
the shocks could resu lt in an afterglow, and t he lum inosity and characteristic photon energy 
increases with time. iKornreich fc Lovelacd (120081 ) discussed that the SMBH displacement 
may give rise to observable non-axisymr aetries in the morphology and dynamics of the stel- 
lar and gaseous disk of the host galaxy. Ide la Fuente Marcos fc de la Fuente MarcosI (120081 ) 
examined the influence of a runaway SMBH passage through intergalactic medium, and indi- 
cated that the SMBH is ab le to ignite star formation efficiently in the wake of its trajectory. 
Komossa fc MerrittI (120081 ) indicated that a recoiling SMBH carries stars, and the electro- 
magnetic flares from the stars that are tidally disrupted by the SMBH would be observable. 



Recently, iBlecha &: Loebl (120081 ) calculated the trajectory of a SMBH ejected in a 
smooth background potential that includes both a stellar bulge and a gaseous disk (see also 



Vicari. Capuzzo-Dolcetta. fc Merrit 



SMBH as a function of time. iFujita 



2007), and estimated the gas accretion rate onto the 



(l2008l . hereafter Paper I) also calculated the trajectory 



of a SMBH ejected in a realistic background potential of a disk galaxy. We calculated the 
accretion rate of gas onto the SMBH from the interstellar medium (ISM) in the galactic disk, 
and estimat ed the X-ray l uminosity based on a model of a radiatively inefficient accretion 
flow (RIAF: lNarayanll2005l ) . We showed that the luminosity of the SMBH can be comparable 
to or even larger thar i those of ultra-luminous X-ray sources (ULXs) observed in galaxie s 



{Lx > 3xl0^^ergs-^: IColbert fc Mushotzkvlll999l : lMakishima et al.ll2000l : lMushotzkvll2004f ) . 
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However, in that study, the effect of dynamical friction was not exphcitly included. If 
dynamical friction is effective, the SMBH gradually settles down to the galactic center. Since 
the accretion rate depends on the density of the surrounding ISM (see equation [H]), we 
expect that the X-ray luminosity increases accordingly. In this paper, we study the long- 
term orbital and luminosity evolution of a recoiling SMBH in a disk galaxy, considering the 
effect of dynamical friction. This paper is organized as follows. Our models and choice of 
parameters are outlined in § [2l The results of calculations are presented in § [3] and discussed 
in § m Finally, in § El we present our conclusions. 



2. Models 



We consider SMBH mergers in a normal disk galaxy for the sake of simplicity, although 
the galaxy interaction and merger supposedly affect the original disk. However, the process 
of the settling of SMBHs between the galaxy merger and the set-in of effective emission of 
gravitational wa ves from the SMBHs, which leads to the merger of the SMBHs, has not b een 



understood (e.g. iBegelman. Blandford. &: Reelll980l : llwasawa. Funato. &: Makind 120061 ). If 



the time-scale of the SMBH merger is long, the galaxy may be significantly relaxed when 
the recoil occurs (IBlecha &: Loebll2008l ). Anyway, in order to follow the galaxy and SMBH 
mergers self-consistently, we would need ultra-high resolution simulations of galaxy mergers 
that can resolve the settling of SMBHs in the core of the merged galaxy. 

The model of a disk galaxy is the same as that in Paper I. The galaxy potential con- 



sists of three com ponents, which are a lMivamoto fc Naeail (19751) disk, Hernquist spheroid 



(IHernquistl Il990l ) . and logarithmic halo (IBinney fc Tremaine 



l2008f ): 



dist 



disk 



R^ + {a + Vz^ + b^y 



halo 



^sphere 



sphere 



r + C 



w + 



(2) 
(3) 



where R (= a/x^ + y"^) and z are cylindrical coordinates aligned with the galactic disk, and 



yB? + z^. The parameters are the ones of the Galaxy. We take Mdisk 



r 

Msphere = 3.4 X 10^" Mq, o = 6.5 kpc, h = 0.26 kpc, c = 0.7 kpc, d 



1.0 X 10^^ Mq, 
13 kpc, and q = 0.9; fhaio 
is determined so that the circula tion velocity for the total potential , $ = $ disk H-'^* sphere +'^'haio, 



is 220 km s ^ at i? = 7 kpc (see lLaw. Johnston, fc Majewskill2005l ). Contrary to Paper I, we 
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include dynamical friction in the equation of motion for the SMBH (jChandrasekharl Il943 



Binney fc Tremaindl2008l ): 



relA 



(4) 



where v = {vx,Vy,Vz) is the velocity of the SMBH, InAj is the Coulomb logarithm, and 
the suffix i refers to disk, sphere, or halo. The relative velocities are defined as i^rei,disk = 
V - f cir.disk, ■frei.sphere = V, and t^rei,haio = V, whcrc t^cir.disk IS the circulatiou velocity of the 
disk, which is given by 



cir, disk 



9$ 



z=0 



Pdisk 



(5) 

(6) 
(7) 



( iBinney fc Tremaindl2008l ). The densities are given by 

'^'MdiskA aR^ + (a + 3Vz2 + b''){a + Vz^ + b'')^ 

47r ; + (a + ^^2 + 52)2]5/2(^2 + ^2)3/2 ' 

-^^splicrc C 

Pspliere 7 ; ; 

ZTT r[c + ry 

( Miyamoto &: Nagai ll97.4lHernqi"fi^ll99nl : iBinney fc Tremaind 120081 ) . We assume that part 
of the disk consists of the ISM; its density is represented by pisM = /isMPdisk- For most 
mod els we studied, we take /j sm = 0.2, which is based on the observations of the Galaxy 



[e.g. 



Binney fc Tremaind l2008l ) . In equation (jlj), the factor I{Xi) is given by 



J(X,) = erf (X,) - 



2X 



-X? 



(9) 



where Xi = Vre\,i/ {V^ai) and ai is the one-dimensio nal velocity dispersion. For the disk 
velocity dispersion, we set Udisk oc Pdisk(a;, ^ = 0) (ILewis fc FreemanI Il989l ) and fix the 
normalization by assumin g that the disk has a Toomre Q-parameter of 1.5 at i? = 7 kpc 
( jVelazquez fc Whitdll999l ). For the spheroid and halo, we use the common I{Xi) and the 
velocity dispersion is 



'-'^sphere— halo 



cir, sphere 



_L.;2 U/2 
' cir, halo/ 



(10) 



where t'rir, sphere and r^rir halo are the circulation velocities of the spheroid and halo, respectively 
jTavlor fc Bab^boOll ). 

Chandrasekhar's formula for the dynamical friction force (the second term in the right 
side of equation |4]) was derived assuming an infinite, homogeneous, and unchanging back- 
ground. It is obviously not true for a SMBH that is kicked out of the galactic center and 
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orbiting in the disk galaxy. However, A^-body simulations have shown that it can be applied 
to various cases if one chooses the Coulomb loga rithm appropriately. For a S MBH ejected 
from the center of a spherically symmetric galaxy, iGualandris fc MerrittI (120081 ) showed that 
2 < InA < 3 is appropriate. Thus, we set InAgphere = InAhaio = 2.5. For a disk galaxy, such 
A^-body simulations have not been performed as far as we know. Instead of a SMBH, the 
evaluation of Chandr asekhar's formul a has b een made for a dwarf galaxy infalling to a more 
massive disk galaxy. iTaylor fc Babull (120011 ) showed that InAdisk = 0.5 is appropriate. The 
small value of In Ajisk probably reflects the small ratio of the disk scale-height and the size 
of a dwarf galaxy, which corresponds to the ratio of maxin aum and miniuni effec tive impact 
parameters of particles that contribute to the friction force (ITaylor &: Babulll200ll ). However, 
since a SMBH is a point source, we first assume that InAdisk = 2.5, which is the same as 
InAgphere a^d InAhaio, and then change the value to see the influence of the uncertainty on 
results. 

The acc retion of the surrounding gas onto an isolated black hole has been studied by sev- 

eral a uthors JFuiita et al.lll998l : lAgol fc Kamionkowskilbooi iMii fc Tota"iiill2005l : [Mapelh. Ferrara fc Rea 
20061 . and references therein). Most of the previous studies focused on stellar mass (~ lOM©) 
or intermediate mass black holes (IMBHs; ~ 10^ Mq). In this study, we consider the accretion 
on a recoiling SMBH. 

The acc retion rate of the ISM onto the SMBH is given by the Bondi-Hoyle accretion 
jBondilll952l ): 



771 = 2.5nG^ 



^bhPism 



11) 



where mBH is the mass of the black hole, and (= 10 km s~^) is the sound velocity of the 
ISM. The X-ray luminosity of the black hole is given by 



(12) 



where rj is the efficiency. Since the accretion rate is relatively smal l for the mass of the 



black hole, the accretion flow would be a RIAF (e.g. Ilchimarul 119771 : iNarayan &: Yil Il994j : 
Abramowicz et al.lll995l : lYuan et al.l 120041 ). In this case, the efficiency follows rj oc m fo r 
Lx < O.lLEdd, where L^dd is the Eddington luminosity (e.g. lKato. Fukue fc Mineshigelll998l ). 
Therefore, we assume that r] = rj^dd f or m > 0.1mF,HH an d 77 = ?7Edd'^/(0.1mEdd) for ^ < 



O.lmEdd, where mEdd = -Z^Edd/(c^'7Edd) ( Mii fc Totani 20051 ). We assume that r^Edd = 0.1. 



We solve equation (jlj) with Mathematica 6.0 using a command NDSolve. The algorism 
of the integration is automatically chosen 0. We have confirmed that the fractional energy 



http: / /support. wolfram.com/matheniatica/mathematics/numerics/ndsolvereferences. en. html 
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error that arises in the integration of an orbit per cycle is < 10 



3. Results 

The SMBH is placed at the center of the galaxy at t = 0. The SMBH is ejected on 
the x-z plane at t = 0. The parameters of our models (mBH, "^O; InAdisk, and /ism) are 
shown in Table [H In this section, we consider models in which the SMBH does neither 
fall into the galactic center too quickly through dynamical friction nor escape from the 
galaxy (models A1-C4 in Table [T]). The mass of the SMBH is 3 x 10^-3 x 10^ M©, which 



Schodel et al. 


2002; 


Ghez et al. 


2005) 



from ^ = 0° to 90°, where 6 
is vq = 500-800 km s'^ 



0° corresponds to the z-axis. The initial velocity of the SMBH 



We stop the calculation if (i) t = 10 Gyr or if (ii) r < 10 pc and v < lkms~^ is satisfied. 
We define the time when the condition (ii) is satisfied as tdf • Figures [TH3] show the trajectories 
of the SMBHs for models Al, B2, and C3, respectively. Although the SMBHs are ejected 
on the x-z plane at t = 0, they are not confined to the plane because of the circulation of 
the galactic disk and dynamical friction. Figures SHB] show the distance from the galactic 
center (r) and the luminosity of the SMBHs {Lx) for models Al, B2, and C3, respectively 
(0 < t < tdf and 9 = 60°). The distance gradually decreases through the dynamical friction. 
The infall of the SMBHs accelerates as t approaches tdf- The luminosity on average increases 
as the SMBHs spiral into the galactic center, where pisM is large. In Table [H we present the 
maximum distance from the center of the galaxy when 6 = 60° (rmax.eo)- We note that the 
maximum radius is not much dependent on 6. 

Figures [7H9] show the evolutions of \z\, frei,disk, and Lx for < t < O.ltdf for models Al, 
B2, and C3, respectively {6 = 60°). The luminosity of the SMBHs {Lx) increases instanta- 
neously, when they pass the galactic disk. The heights of the spikes in Figure [7b, [S)d, andlHb 
are uneven. This is because the luminosity Lx depends on both frei,disk and pisM (see equa- 
tion [H]), and the latter strongly depends on z. In Tabled we present the maximum X-ray 
luminosity of the SMBHs for 9 = 60° and t < 0.2tdf (-^^max,6o)- The luminosity Lmax.eo is larger 
for larger and smaller vq. In some of the models of > 1 x 10^ Mq, Lmax,60 reaches 
3 X lO^^ergs"^, which is the Eddington luminosity of a stellar mass b l ack hole (~ 20 M(^) and 



is often used as a threshold of ULXs (IColbert fc MushotzkyI Il999t iMakishima et al.l 12000 



MushotzkyI 120041 ). When mnn and fo are given and 6 is not fixed, the X-ray luminosity 
before the SMBHs are affected by dynamical friction tends to be larger when 6 is closer to 
90°, because their trajectories are included in the galactic disk, where pisM is large. However, 
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the tendency is not clear when 9 < 80°, because their orbits are scattered in the asymmetric 
potential of the galaxy. Thus, the maximum luminosity does not much depend on 9. Fig- 
ures [7HH] indicate that when the SMBH is especially bright, the relative velocity between the 
SMBH and the surrounding ISM (or stars) is frei,disk ~ Vc\r{^ 220kms~^), which means that 
the SMBH passes the apocenter of its orbit [v ~ 0) close to the galactic plane. It could be 
used as a clue to find the traveling SMBH observationally, if atomic line emission associated 
with the X-ray source is detected and the velocity is estimated through the Doppler shift. 

In Table m we present the average of tdf, which is referred to as (tdf); we calculate 30 
orbits and corresponding tdf by changing 9 from 3° to 90° by 3° at a time, and average t^i 
by 6', weighting with sin^. Table [1] shows that (tdf) ^ 0.1 Gyr for models A1-C4, and that 
(tdf) is smaller for larger rriBH and smaller Vq. 

Following Paper I, we estimate the probability of observing SMBHs with luminosities 
larger than a threshold luminosity Lth, assuming that SMBHs are ejected in random direc- 
tions at the centers of galaxies. For each model, we calculate 30 evolutions of the luminosity 
by changing 9 from 3° to 90° by 3° at a time. Then, we obtain the period during which the 
relation Lx > i^th is satisfied for each 9, and divide the period by tdf- This is the fraction 
of the period during which the black hole luminosity becomes larger than Lth- We refer to 
this fraction as f{9). We average f{9) by 9, weighting with sin^, and obtain the probability 
of observing SMBHs with Lx > Ivth- In Table [H we present the probability P^em when 
Lth = 3 X 10^^ erg s^^ For models A1-C4, Pgesg = 0.0018-0.58. 

We also estimate the age-corrected probability of observing SMBHs with Lx > Lth = 
3 X 10^^ erg s~\ which is obtained by averaging min[tdf, iage]/(6')/iagc by 9, weighting with 
sin 6*, where tage is the age of a galaxy and we assume that tage = 10 Gyr. We refer to the 
age-corrected probability as -Pjesg and show it in Table [TJ 

4. Discussion 

We have found that a SMBH that had been ejected from the center of a disk galaxy 
could be observed in the galactic disk with an X-ray luminosity of Lx ^ 3 x 10^^ erg s~^. 
The luminosity gradually increases as the SMBH settles down to the galactic center through 
dynamical friction. 

In § [3|, we follow the evolution until the SMBH spirals down to r = 10 pc. However, 
if r is too small, the SMBH cannot be discriminated from the one that would have been 
sitting at the galactic center without being affected by a recoil. Therefore, we estimate the 
probability of observing SMBHs with Lx > Lth = 3 x 10^^ erg s~^ and r > 1 kpc, and call 
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it Pi,3e39. We also calculate the time age-corrected one (-Pi,3e39). Since Pi,3e39 and Pi,3e39 are 
derived by adding another condition r > 1 kpc to P3e39 and -P3e39 5 respectively, it is natural 
that Pi,3c39 < Ae39 and Pi,3e39 < ^3c39 (Table [I]). 

As is mentioned in § [21 dynamical friction of a massive point particle orbiting in a 
disk galaxy has not been studied very much. Thus, there is some uncertainty about the 
Coulomb logarithm we should take. Therefore, we change the value of InAdisk to estimate 
the uncertainty. Models C'2 and C'3 are respectively the same as models C2 and C3 except 
for InAdisk- For these models, we set InA^isk = 1-5. Table [T] shows that there is not much 
difference between the results of models C2 and those of C'2. This is because the maximum 
distances to the apocenters are < 1 kpc, where the spheroidal component is dominant, and 
the SMBH is not much affected by the dynamical friction from the galactic disk. On the 
other hand, P3e39, -P3e395 -Pi.sesg, and A.sesg for models C3 and C'3 are significantly different, 
because the SMBH is ejected outside the spheroid. The differences are especially made by 
that of the orbits of 6* ~ 90°. When 9 ~ 90°, the SMBH is ejected in the galactic disk. If 
the dynamical friction from the disk is very effective, the SMBH moves along with the disk 
(■^rei.disk ~ 0) and docs not easily fall into the galactic center. This actually happens for 
model C3 (tdf > 10 Gyr when 9 = 90°; Figure [TUk) . In this case, the SMBH continues to 
accrete the ISM in the disk and is bright for a long time (Figure ITOb). Since this SMBH is 
very bright {Lx ^ 10^^ erg s~^), it could be easily observed if such SMBHs actually exist. 
For model C'3, the dynamical friction from the disk is not strong enough to hold back the 
SMBH from the infall even when 9 = 90°. 

We also consider the uncertainty of the ISM fraction /ism- Model C"3 is the same as 
model C3 but for /ism = 0.1. For parameters we adopted, the accretion efficiency is m < 
0.1 mEdd in most cases. Therefore, we obtain Lx oc rjm oc oc pjgj^ oc /j^j^ (equations [11] 
and [12]), which means that the X-ray luminosity in model C"3 is one fourth of that in 
model C3. Accordingly, P3039, -P3e39, Pi,3e39, and Pi,3e39 in model C"3 are smaller than those 
in model C3, respectively. However, the differences are not large, because Lx changes rapidly. 

For comparison, we also investigate a model with a smaller initial velocity (model bO), 
because we consider a Milky- Way type galaxy, which is generally expected to experience mi- 
nor mergers rather than major mergers. In such cases, large recoil velocities as adopted above 
would not be common. A model with a larger SMBH mass is also considered (model d3). 
Figures [11] and [T2] show the trajectories of the SMBHs for models bO and d3, respectively. 
Their ejection angles are 9 = 60°. In these models, Vq is too small (model bO), or 
is too large (model d3) for the SMBH to be ejected from the spheroidal component of the 
galaxy. Thus, it would be difficult to recognize them as recoiling SMBHs, if their host galax- 
ies are moderately distant. The SMBHs are almost confined to the x-z plane, because the 
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dynamical friction from the spheroidal component overwhelms that from the galactic disk. 
Table [1] shows that Psesg for models bO and d3 is relatively large because of small fo and 
large rriBK, respectively. The SMBHs set back to the galactic center in only several orbital 
periods (tdf ~ 0.01 and 0.07 Gyr, respectively). 

Since we have included the effect of dynamical friction when we consider the evolution 
of Lx, we can constrain the probability to find SMBHs with Lx > -^^th more precisely than 
Paper I. It has been estimated that for comparable mass binaries with dimensionless spin 
values of 0.9 , only ~ 10% of all mergers are e xpected to result i n an ejection speed of ~ 500- 



800 km s ^ ( ISchnittman fc Buonannd 120071 : iBaker et al.ll2008l ). Since the ejection speed is 



smaller for mergers with large mass ratios and smaller spin values, the actual fraction would 
be smaller. Although we consider the mergers of black holes with the masses currently 



Enoki et al. 


2004; 


Micic et al. 


2007) 



The number of such mergers that a galaxy has undergone would be < 1. Thus, since 
PseSQ, -Pi,3e39 ^0.1 (Table [T]), the probability that a disk galaxy has a traveling SMBH with 
a luminosity comparable to or larger than that of ULXs is < 1 x 10~^. 

Since the probability is not so large, extensive surveys would be required to find the 
SMBHs running in the galactic disks. In the future, it would be interesting to study whether 
the probability is larger than that of finding SMBHs immediately after the ejection from the 



galactic centers with velocities of > 1000 km s~^(e.g. iLoebl 120071 ). Since the SMBHs ejected 
into galactic disks are very bright (Figure [TUb). they could be observed even in distant 
galaxies. As was discussed in Paper I, observations in bands other than X-rays would also 
be useful to detect the SMBHs orbiting in disk galaxies and discriminate them from IMBHs. 



5. Conclusion 

We have investigated the trajectory of a SMBH ejected from the galactic center through 
the emission of gravitational waves at the merger of two black holes. We included the effect 
of dynamical friction. For a disk galaxy comparable to the Galaxy, the orbit decays on 
a time-scale of > 10^ yr if the initial velocity of the SMBH is ~ 500-800 km s~^ and the 
mass is ~ 10'' Mq. The SMBH accretes the surrounding ISM when it passes the galactic 
disk. Since the accretion rate is larger when the relative velocity between the SMBH and 
the ISM is smaller, the accretion rate is the largest when the SMBH passes the apocenter 
of its orbit that reside in the galactic disk. Assuming that the accretion flow is a RIAF, we 
estimated the X-ray luminosity of the SMBH. We found that the X-ray luminosity can reach 
Lx ^ 3 X lO^^ergs"^, which is comparable to or even larger than those of ULXs. In particular. 
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the X-ray luminosity would reach Lx ^ lO^^ergs"^, if the SMBH is ejected into the galactic 
plane. Since the probability of finding the traveling SMBHs with Lx ^ 3 x 10^^ erg s~^ in a 
disk galaxy is < 0.01, extensive surveys would be required to find them. 
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Sports, Science, and Technology of Japan (20540269). 
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2. — Same as Figure [1] but for model B2. 
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Fig. 3. — Same as Figured] but for model C3. 
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Fig. 4. — (a) The distance from the galactic center and (b) the luminosity of the SMBH for 
model Al for < t < tdf when 6 = 60°. 
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Fig. 5. — Same as Figure H] but for model B2. 




Fig. 6. — Same as Figure H] but for model C3. 
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Fig. 7. — (a) The distance from the galactic plane, (b) the relative velocity between the 
SMBH and the galactic disk, and (c) the luminosity of the SMBH for model Al for < i < 
0.1 idf when ^ = 60°. 
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Fig. 8. — Same as Figure [7] but for model B2. 
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Fig. 9. — Same as Figure [7] but for model C3. 
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Fig. 10. — (a) The trajectory in the galactic disk plane {z = 0) and (b) the luminosity of 
the SMBH for model C3 for < i < 10 Gyr when = 90°. 
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Fig. 11. — The trajectory of the SMBH for model bO for < t < t^f when 6 = 60°. 
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Fig. 12. — Same as Figure [TT] but for model d3. 




